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Abstract 

Presented is a mathematical model, derived from the 
Navier-Stokes equations of momentum and continuity, 
which may be accurately used to predict the behav- 
ior of conventionally mounted pneumatic sensing sys- 
tems subject to arbitrary pressure inputs. Numerical 
techniques for solving the general model are developed. 
Both step and frequency response lab tests were per- 
formed. These data are compared against solutions of 
the mathematical model. The comparisons show ex- 
cellent agreement. The procedures used to obtain the 
lab data are described. In-flight step and frequency re- 
sponse data were obtained. Comparisons with numer- 
ical solutions of the mathematical model show good 
agreement. Procedures used to obtain the flight data 
are described. Difficulties encountered with obtaining 
the flight data are discussed. 

Nomenclature 

cross-sectional area of pneumatic tubing 
sonic velocity 
polytropy coefficient 
pressure tubing cross-sectional diameter 
spatial index for discretization grid 
temporal index for discretization grid 
length of pneumatic tubing 
maximum number of temporal 
gridpoints 

maximum number of spatial gridpoints 
pressure function 
ambient pressure 
i ih spatial, k th temporal pressure 
gridpoint 

inverse pressure state, 
acoustic resistance of pneumatic 
configuration 

Reynolds number of tubing flow 
time coordinate 
velocity function 
i th spatial, k th temporal velocity 
gridpoint 
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enclosed transducer volume 
spatial coordinate 
ratio of specific heats for air 
discretization parameter, (^f-) 2 
temporal grid stepsize 
spatial grid stepsize 
discretization parameter, 1 + RAt/p 
dynamic viscosity of air 
polytropy heat transfer parameter 
density function 


Introduction 


Recent advances in aircraft performance and ma- 
neuver capability have dramatically complicated the 
problem of flight control augmentation. Some control 
system designs now require that aerodynamic parame- 
ters such as angle of attack, dynamic pressure, or wing 
pressure loads be fed back to the control system. The 
use of aerodynamic parameters as control system feed- 
backs requires the measurements be recorded with ac- 
curacy and high fidelity. Since most aerodynamic pa- 
rameters must be pneumatically sensed, this is a diffi- 
cult problem. The primary difficulty in obtaining high- 
frequency pressure measurements is pressure loss ow- 
ing to frictional attenuation within the sensing system. 
Typically, most of the frictional loss comes within the 
pneumatic tubing used to transmit pressure impulses 
from the surface of the aircraft to the measurement 
transducer. To avoid pneumatic attenuation, designers 
have sought to mount the pressure sensor at the surface 
of the aircraft (in situ mounting) or, when this is im- 
possible, at least restrict the tubing length to just a few 
inches. In some cases this is a viable solution. How- 
ever, when many pressures must be measured in a small 
surface area, as with a hemispherical airdata sensor, 
one simply cannot crowd enough pressure transduc- 
ers into the available space. As a result, the designer 
must accept whatever compromise is available. De- 
pending upon the aircraft, this means running sizeable 
lengths at pneumatic tubing from the surface to the 
pressure transducer. 

The emphasis of this paper is on the development of 
a general numerical technique for accurately predicting 
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pneumatic attenuation errors when the pressure trans- 
ducers cannot be mounted in situ. Accurate prediction 
of pneumatic attenuation errors allow control system 
designers to provide for adequate feedback-loop robust- 
ness. The model is formulated to allow for arbitrary 
pressures to be input to the sensing system. Numeri- 
cal techniques for solving the mathematical model are 
developed. The model is verified by comparing both 
laboratory and flight data. 

Background 

As mentioned in the Introduction, for many flight 
test applications pressure-sensing devices cannot be 
easily or practically mounted flush to the aircraft skin 
(in situ mounting). To transmit pressure change at 
the surface to the transducer, a length of connective 
tubing is used. Pressure variations at the surface prop- 
agate as waves from the upstream end through the con- 
nective tubing to the transducer. The wave propaga- 
tion is damped by frictional attenuation along the walls 
of the tubing. The wave damping manifests itself as 
spectral attenuation of the pressure response and pro- 
duces both a magnitude attenuation and a phase lag. 
When the wave reaches the downstream end of the tub- 
ing, it is reflected back up the tube and may either 
damp or amplify incoming pressure waves leading to 
spectral distortion. 

A considerable body of information concerning the 
effects of pneumatic attenuation is available. Early 
attempts made use of acoustical-electrical analogs to 
derive time-dependent low-order linear models to ap- 
proximate pneumatic attenuation. 1,2 Such analyses, al- 
though of some predictive value, are not based on rig- 
orous aerodynamic principles and have little generality. 
Later analyses were based on more rigorous approaches 
using the equations of momentum and continuity. Such 
analyses lead to response models for specialized pres- 
sure inputs, such as step inputs, 3 ramp inputs, 4 and 
sinusoidal inputs. 5,6 These analyses, while accurately 
predicting the pressure response for their prescribed 
types of inputs, are unable to predict the behavior 
of pneumatic systems subjected to arbitrary pressure 
inputs. A more general mathematical model must be 
developed. 


Idealized Configuration Geometry 

The sensor configuration is modeled as a straight 
cylindrical tube with an axisymmetric volume attached 
to its downstream end. The tube represents the trans- 
mission line from the surface to the sensor, and the 
attached volume represents the internal volume of the 
pressure transducer. The tube is considered to be of 
constant radius R , with length L. The attached volume 
V is assumed to be a constant. A longitudinal coordi- 
nate x is measured from the upstream end of the tube, 
and a time coordinate t is measured forward from some 
initial time J 0 - Pressure, density, temperature, and the 
radial average flow velocity within the tube are consid- 
ered to be functions of only x and t. These quantities 
are represented by the symbols P(x , t ) , p(x, t ) , T(x, <), 
and U(x 1 t) ) respectively. The idealized configuration 
is depicted in Fig. 1. 

Navier-Stokes Equations 

The Navier-Stokes equations, contingent on the 
aforementioned assumptions, are given by Stephens 
and Bate 7 as 

(i) continuity 

| 01 _ 0 /jx 

dt dx 


(ii) momentum 


dU(x,t) 
P dt 


+ 


dP(x,t) 

dx 


+ RU (x, t) = 0 


( 2 ) 


The term R in equation (2) is referred to as the 
acoustic resistance and represents the damping effects 
of viscosity. For laminar flow the value of R may be 
shown to be 



(3) 


where /i is the local dynamic viscosity of air. For tur- 
bulent flow conditions, the friction coefficient may be 
estimated using Blasius’ empirical law 8 to give 


= 0.1582/2# 


D 2 


(4) 


Derivation of Mathematical Model 

This section presents the idealized configuration to 
be analyzed first. Next the Navier-Stokes equations of 
momentum and continuity are presented. The equa- 
tions of momentum and continuity are coupled to 
give a single nonlinear partial differential equation for 
pressure in terms of position and time. The nonlin- 
ear equation is linearized to give a damped wave equa- 
tion. Finally, the initial and boundary conditions 
are developed. 


Empirical data presented in Schlicting 8 indicate that 
pipe flow transitions from laminar to turbulent for 
Reynolds numbers above 2300. 

Reduction to Nonlinear Wave Equation 

If the density behind the pressure wave front at any 
given time is assumed to be independent of position, 
and the wave process is assumed to be polytropic 4 
where 

£ = c ° <») 
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equations (1) and (2) may be combined to form a single 

nonlinear pressure wave equation P(0,t) = Pi(t) (11) 


d 2 P(x,t) R dP(x,t) 2 d 2 P{x,t) 
di 2 + p(t) dt y C dx 2 
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+ 4P{*,t) 


’ dP{t,t) ] 2 

3t 


( 6 ) 


where 7 is the ratio of specific heats, and c is the local 
sonic velocity. For the current analysis, it is assumed 
that the flow velocities are low, hence c is approxi- 
mately constant. The parameter £, known as the heat 
transfer parameter, satisfies 1 < £ < 7. If air is the 
fluid medium, and pressure changes are not large, then 
£ = 7 = 1.4 is sufficiently accurate. 

Linearization of Nonlinear Wave Equation 

Equation (6) may be effectively linearized by noting 
that 

P(x,t) = C —p(x,t) 

Substituting into the last term on the right-hand side 
of equation (6) and regrouping gives 


d 2 P(x,t) 

dt 2 
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dPjxJ) = j a fl 2 P(«,«) 
dt dx 2 

(7) 

Assuming that input pressure rates are not large, an 
order of magnitude analysis for a typical pressure sens- 
ing system 4 indicates 
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Thus equation (7) is effectively linearized to give 
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Equation (9) is simply the classical wave equation 
with time variable parameters. The pressure variation 
within the tube can be visualized as the propagation 
of a longitudinal compression wave. 

Development of Initial and Boundary 
Conditions 

To complete the mathematical model, the initial and 
boundary conditions must still be defined. The initial 
conditions will be defined by assuming that at time fo, 
the system is at rest. As a result the initial pressure 
and velocity at all stations are prescribed to be 


P(x,0) = P o , U(x, 0) = 0 (10) 


The boundary condition at the upstream of the tube 
is prescribed 


The boundary condition at the downstream end of the 
tube may be obtained by satisfying the equations of 
continuity and momentum at the downstream end to 
give 

= w) 

P L (t) V 

and, 


( 12 ) 


P^P- + RU L (t) + 


dP(x, t) 


dx 


0 


(13) 


where A c is the cross-sectional area of the tube. 

Numerical Solution of Pneumatic 
Attenuation Model With 
Arbitrary Inputs 


The mathematical model given by equations (9), 
(10), (11), (12), and (13) is to be solved numerically 
using an implicit-differencing technique. The pres- 
sure function is approximated by a series of spatially 
and temporally discrete gridpoints, with partial deriva- 
tive operators being approximated by finite differences. 
For each temporal recursion, the entire spatial grid 
is solved, and the differential equation for the down- 
stream boundary is updated by numerical integration. 

Finite Difference Approximation of Wave 
Equation 

As previously mentioned, the pressure function will 
be approximated by a series of discrete gridpoints 


P(x,t) = p ik 


and 


U(x,t) = u ik 


with partial derivatives being approximated by differ- 
ence operators; that is, first derivatives are approxi- 
mated by 

dP(x,t) _ Pit+i - Pi „ 
dt At 

and 

dP(x } ^) Pi + l k+ i ~ Pik + I 


dx 


Ax 


(14) 

(15) 


The index i represents the i ' th spatial gridpoint, and 
the index k represents the Ar'th temporal gridpoint; 
while Ax and At represent the spatial and temporal 
distances between the discrete gridpoints. The differ- 
ence operators are substituted into equation (9) and 
rearranged to give the recursive relationship 


— ^Pt+i*+i+(^+ 2 A)p , 4+1 — Ap ,_ lk+1 = ~(l+%, 

(16) 
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where A = (cAt/Ax) 2 and 6 = 1 4- [/2Af //>(*)]. Equa- 
tion (16) must be satisfied for 

s — 1, 2, ... M, M — ,and 
Ax 

*=1,2, ... AT = Lli® 


At temporal gridpoint k+1, the spatial grid may be 
written in matrix form as the symmetric tridiagonal 
system 


' 6 + 2\ 

-A 

o 

-A 

6 + 2A 

-A 

0 

-A 

6 + 2A 

0 

0 

o 

0 

0 

o 

0 

0 

o 


0 0 O' 


P H+i 

0 0 0 


P 2*+i 

0 0 0 


P3j.+i 

• • • 6 + 2A -A 0 


PM- 3*+, 

-A 6 + 2A -A 


PM-2 k+l 

0 -A 6 + 2A 


. PM- 1* +1 . 




Pl*-1 


•*Po»+, 

(1 + £)P2» 


P 2.-, 


0 

(1 + ^)P3* 


P3*_, 

+ 

0 

(1 +6)pM-3 k 


PM- 3*_, 


0 

(1 + S)PM-2 k 


PAf — 2»_i 


0 

. (1 + 6 ) pm - i „ . 


. PM-U-i 


. *PN k+l . 


Evaluation of Initial and Boundary Conditions 

The initial condition and the upstream boundary 
condition are given explicitly by equations (10) and 
(11). The downstream boundary condition, given by 
equations (12) and (13), must be evaluated by integrat- 
ing with respect to time. Unfortunately, the boundary 
condition given by equation (12) is unstable for positive 
values of £/l(<). This instability requires reformulat- 
ing the downstream boundary condition in terms of an 
inverse pressure state. This reformulation, which 
will circumvent instability problems, follows. Solving 
equation (13) for Ul(1) and substituting into equa- 
tion (12) gives 


and substituting into equation (18). Regrouping and 
integrating from temporal indices k to k + 1 gives 


Pk+i 




[ 1 + 1 

£A c At r 7Ar(u Nfc+1 -t4^ Jc ) ( 1 1 

kvtzzI- " c *lt + 1 J 


+ 


[ 1 + RV^fgPN- U+i 1 


(19) 


Equation (19) is a stable reformulation of equation 
(12) and may be used to evaluate the downstream pres- 
sure at each temporal gridpoint. 

An implicit algorithm is used for solving equa- 
tion (9), subject to the constraint of equations (10), 
(11), (13), and (19). The method requires iteration to 
convergence at each temporal gridpoint. Each tempo- 
ral recursion consists of three separate steps: a predic- 
tion step, a correction step, and an iteration step. The 
details, since they are mostly concerned with bookkeep- 
ing and stability checking, will not be presented. 


Verification of Mathematical Model 


At this point, it is appropriate to switch emphasis to- 
wards verification of the model. Empirical verification 
of the mathematical model is presented in two steps in 
this section. First, numerical solutions are compared 
to laboratory data. Specifically, numerical solutions 
are compared to data obtained from step response and 
frequency response laboratory experiments. Second, 
numerical solutions are compared to actual flight data 
obtained from a specifically designed experiment. The 
flight data comparisons aid in establishing the appli- 
cability of the mathematical model to realistic pressure 
measurement configurations. 

Laboratory Step Response Tests 

Step Response Test Equipment and Proce- 
dures. This test assembly was designed to demon- 
strate the effects of pneumatic attenuation upon step 
inputs. The general procedure consisted of comparing 
the response of an in situ mounted reference transducer 
to the response of a test transducer mounted at the end 
of a sizable length of pneumatic tubing. A ^ ^4 step 
was input to the system using an evacuated chamber 
and an electrically actuated solenoid valve. In order to 
simulate altitude effects, the experiment was designed 
to allow the tests to be performed at a variety of am- 
bient pressures. 


PL(t) VR [ P dt dx 


(18) 


Equation (18) may be transformed to a linear equation 
by defining an inverse pressure state 


p-\t) 


i 

PlW 


Results of Step Response Tests. The results 
from the step response tests are now compared against 
numerical solutions of the mathematical model. Fig- 
ures 2 and 3 present sample step response time his- 
tories. In these figures, the solid line represents the 
response of the test transducer as measured in the lab, 
the short dashed line represents the numerically com- 
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puted response of the test sensor, and the hashed line 
represents the step input. Figure 2 depicts the response 
of a test configuration with 40 ft of 0.081 in. steel tub- 
ing at 2300 ft simulated altitude. Figure 3 depicts the 
response of a test configuration with 20 ft of 0.081 in. 
steel tubing at 20,000 ft simulated altitude. The agree- 
ments are excellent. 

Laboratory Frequency Response Tests 

Frequency Response Test Equipment and 
Procedures. This test assembly was designed to 
demonstrate the effects of pneumatic attenuation upon 
sinusoidal inputs. The general procedure consisted of 
comparing the relative frequency responses of an in situ 
mounted transducer to a transducer mounted at the 
end of a test section of pneumatic tubing. A piston- 
actuated, electromechanical acoustical amplifier was 
used to input sinusoidal waves of prescribed frequen- 
cies and amplitudes to the measurement configuration. 
The base pressure within the amplifier cylinder cavity 
could be raised or lowered to simulate a desired alti- 
tude. Tests were repeated for each of the tubing test 
sections at a variety of base pressures. 

Results of the Lab Frequency Response Tests. 
Selected results of the frequency response tests are now 
compared against numerical solutions of the mathe- 
matical model. Figures 4 and 5 present sample Bode 
plots resulting from the lab tests. Both magnitude 
and phase angle plots are presented. In these figures, 
the solid line represents the frequency response of the 
test transducer as measured in the lab, and the short 
dashed line represents the numerically computed fre- 
quency response. Figure 4 depicts the frequency re- 
sponse of a test configuration with 24.5 in. of 0.046 in. 
steel tubing at 2300 ft simulated altitude. Figure 5 
depicts the frequency response of a test configuration 
with 24.5 in. of 0.046 in. steel tubing at 20,000 ft sim- 
ulated altitude. As with the step response data, the 
matches are excellent. 

The comparisons of this section indicate that the 
ability of the model to predict the response of sim- 
ple pressure measurement configurations is quite good. 
The model is equally accurate for step and frequency 
response predictions. 

Comparisons of Flight Data to Numerical 
Solutions of Mathematical Model 

With the accuracy of the mathematical model in 
predicting the behavior of simple configurations un- 
der a controlled laboratory environment being demon- 
strated, it now becomes pertinent to demonstrate the 
applicability of the model to an actual flight test sens- 
ing system. Flight tests which approximated both step 
and frequency response experiments were performed. 

Flight Test Equipment and Procedures. The 
following section describes the configuration that was 


used to obtain the flight test data. Basically, the con- 
figuration consisted of a matrix of five pressure sen- 
sors with ports clustered at the 10 percent chord and 
20 percent span on the upper surface of the right-hand 
wing of an F-15 aircraft. The flight test configuration, 
located on the right-hand wing, is shown in Fig. 6. 

Two of the ports were configured to serve as refer- 
ence sensors. The other three ports were configured to 
serve as test sensors. The reference transducers were 
mounted so as to produce little pneumatic attenuation 
in the pressure measurements. The test sensors were 
mounted at the end of tubing sections similar to those 
used in the lab tests. The tubing sections were in- 
terchangeable to allow for the testing of various tub- 
ing sections. 

Using the following technique, step response data 
were obtained at a variety of altitudes from 10,000 to 
40,000 ft. First, the aircraft was flown at a high sub- 
sonic speed until a shock formed downstream of the test 
matrix. At this point the pilot would perform a windup 
turn (to the left), thereby increasing the local velocity 
on the surface of the right-hand wing. This caused the 
shock wave to move forward until eventually it passed 
over the sensor matrix. The effect was to introduce a 
sharp step-like increase in the local pressure. As the 
pilot rolled out of the turn, the local velocity on the 
right-hand wing dropped and the shock retreated. The 
effect was to introduce a sharp step-like drop in the 
local pressure. 

Although pure frequency response data were impos- 
sible to obtain in-flight, the following technique proved 
useful in obtaining a broadband fairly white input to 
the test matrix. At moderate subsonic speeds, the pi- 
lot would roll into a windup turn to the right. At this 
point the pilot pulled the aircraft into a high angle of 
attack. The combination of lateral g’s and high angle 
of attack caused mild leading-edge separation. Noisy 
vortices, which would subsequently impinge upon the 
sensor matrix, were shed. The spectral content of these 
vortices was broadband and (beyond 20 Hz) nearly 
white. Because the vortices were random, no phase 
coherence could be established between the reference 
and test sensors. 

Results of the Flight Tests. Selected results of 
the flight tests are now compared against numerical 
solutions of the mathematical model. Step response 
data are presented first. The frequency response data 
will then follow. 

Sample step response time histories are presented in 
Figs. 7 and 8. In these figures the solid line repre- 
sents the measured response of the reference sensor, 
the short dashed line represents the numerically com- 
puted response of the test sensor, and the hashed line 
represents the measured response of the test sensor. 
Depicted are the responses of 2-ft sections of 0.02 and 
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0.04 in. diameter steel pnuematic tubing at an altitude 
of 39,800 ft. Again, the agreements are outstanding. 

Selected results of the frequency response tests are 
now compared against numerical solutions of the math- 
ematical model. Figures 9 and 10 present sample Bode 
plots in which the relative magnitudes of the test and 
reference sensor measurements are compared against 
the numerically computed frequency response of the 
test sensors. In these figures the solid line represents 
the frequency response of the test transducer (rela- 
tive to the reference transducer) as measured in-flight, 
and the hashed line represents the numerically com- 
puted frequency response of the test sensor. Figure 9 
shows the frequency response of a test configuration 
with 24.5 in. of 0.04 in. steel tubing at 40,000 ft al- 
titude. Figure 10 depicts the frequency response of a 
test configuration with 24.5 in. of 0.02 in. steel tubing 
at 40,000 ft altitude. The comparisons are good. 

Conclusions 

The mathematical model, derived from the Navier- 
Stokes equations expressed in one spatial dimension, 
is essentially a damped wave model. Due to the as- 
sumptions required in linearizing the model, it is valid 
for predicting the behavior of pneumatic measurement 
configurations in which the rate of change of the pres- 
sure within the configuration is not large. Since this 
is the case for most pneumatic measurement systems, 
the model has general applicability. 

Comparisons of the mathematical model to both lab 
and flight step and frequency response data indicate 
that its predictive ability is excellent. The model is 
capable of predicting the response of a simple mea- 
surement configuration for any sort of input. It is thus 
directly applicable to flight data. 

By applying the boundary conditions in the reverse 
direction (letting the downstream pressure be directly 
measured and solving continuity and momentum at the 
upstream end), it is possible to invert the mathematical 
model to give a compensation routine that may be used 
to adjust for the effects of pneumatic attenuation. The 
compensation technique still has some problems with 
numerical stability and is being further developed at 
this time. 

There are many difficulties with obtaining in situ 
pressure measurements under a flight environment. 


Primary difficulties lie with controlling the environ- 
ment of the sensor, ensuring survivability, and damp- 
ing out structural noise on the sensor. At this time no 
adequate technology exists which allows one to prac- 
tically and cheaply obtain high-quality pressure mea- 
surements using in situ mounted sensors. Convention- 
ally mounted sensors with short lengths of pneumatic 
tubing have been shown herein to give excellent results 
and are not subject to many of the difficulties encoun- 
tered with in situ mounted sensors. Using the devel- 
oped code, the effects of pneumatic attenuation can be 
identified. Having identified the effects of pneumatic 
attenuation, the control designer can provide for suf- 
ficient robustness in loops that use the pneumatically 
derived data as feedbacks. 
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Fig . 1 Schematic of pressure sensing de- 
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Fig. 2 Step response data: comparisons of 
lab to mathematical model for L = 40 ft, 

D = 0.081 in., and altitude ~ 2300 ft. 
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Fig. 3 Step response data: comparisons of 
lab to mathematical model for L = 40 ft, 

D = 0.081 in., and altitude — 20,000 ft. 
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Fig . 4 Frequency response data: compar - 

isons of lab to mathematical model for 
L = 24-5 in., D = O.O 46 in., and alti- 
tude = 2800 ft. 
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Fig. 5 Frequency response data: compar- 
isons of lab to mathematical model for 
L = 24.5 in., D = O.O 46 in., and alti- 
tude = 20,000 ft. 


0.04-in. x 
diameter \ dla ™ et ® r 
porty y° rt 






reference sensor 



To aircraft 
T PCM system 


I ^Absolute 

pressure gage 


To pulse code 
modulation 
system 


8130 


Fig. 6 Experimental configuration used for flight tests. 
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Fig. 7 Step response data: comparisons of 
flight results to mathematical model for 
L = 2 ft, D = 0.020 in ., and altitude = 
39 , 800 ft. 



Fig . 8 Step response data: comparisons 

of flight results to mathematical model for 
L = 2 ft, D = 0.040 in., and altitude = 
39,800 ft. 
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Measured flight data 



Fig . 9 Frequency response data: compar- 

isons of flight data to mathematical model 
for L = 24-5 in., D = 0 . 0 4 in., and alti- 
tude = 40,000 ft. 
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Fig. 10 Frequency response data: compar- 
isons of flight data to mathematical model 
for L = 24-5 in., D = 0.02 in., and alti- 
tude = 40,000 ft. 
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